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Longitudinal optical phonons have been used to interpret the electronic energy relax- 
ation in quantum dots and at the same time they served as a reservoir, with which the 
electronic subsystem is in contact. Such a phonon subsystem is expected to be passive, 
namely, in a long-time limit the whole system should be able to achieve such a stationary 
state, in which statistical distributions of both subsystems do not change in time. We pay 
attention to this property of the LO phonon bath. We show the passivity property of the so 
far used approximations to electronic transport in quantum dots. Also we show a way how 
to improve the passivity of LO phonon bath using canonical Lang-Firsov transformation. 
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1 Introduction 

Semiconductor nanostructures, like quantum dots, appear to bring new proper- 
ties due to their restricted size in all three dimensions. The size of them is often not 
much larger that the charge carrier's mean free path. Electronic energy relaxation 
in quantum dots does not seem a fully resolved question, both in experiment and 
theory pp. The rapidity of electronic relaxation may be ascribed to the coupling 
of electrons with longitudinal optical (LO) phonons 1211^, although the relaxation 
effect in quantum dots was also interpreted with help of Auger mechanism. Besides 
this, theory of polarons in quantum dots has developed recently, which concept, 
together with the finite lifetime of LO phonons, has became another way of expla- 
nation of the rapid relaxation of electronic energy in quantum dots 4 . 

The quantum transport theory, building in the concept of nonequilibrium Gre- 
en's functions 0, or alternatively on using the theory of nonequilibrium statistical 
operator [^, used the self-consistent Born approximation to the electronic self- 
energy. This approximation to the quantum transport of charge carriers in quantum 
dots was able to provide explanation of the rapidity of electronic relaxation |21 E] , 
and of some properties of the optical transitions, like the line-shape [7], but it also 
recently provided arguments in favor of electronic level occupation up-conversion 
and incomplete depopulation in quantum dots, which effect appears in remarkable 
agreement with a number of experimental papers (see e.g. [51^]). The self-consistent 
Born (SCB) approximation to electronic self-energy, which we are going to use 
here, means that the electron is bound to a cloud of virtual LO phonons. It is 
already at this level that we consider explicitly the motion of the reservoir and 
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its entanglement with the electronic subsystem. This optical phonon cloud can 
be viewed as coherent phonons coupled to an electron. In analogy with coherent 
light, this coherent phonons may act as an effective classical time-dependent field, 
in which the electron moves. Due to this feature we can meet with the off-shell 
property of the collision integral of the electronic quantum transport equation in 
the SCB approximation. 

Because this self-consistent Born approximation to the electron quantum trans- 
port equation leads to an off-shell formulation of the transport equation, a question 
naturally arises concerning an overall energetic balance in the processes of electronic 
energy relaxation via electron-LO-phonon scattering in quantum dots. Within this 
subject we can raise questions about the passivity of the phonon reservoir. In par- 
ticular, we may ask, whether the whole system of electrons and phonons, can evolve 
to a long-time limit, such that in this limit the density matrices of both subsystems 
are in a steady state. If this happens to be the case, we can speak about a passivity 
of the phonon reservoir. Naturally, this passivity property is a desired character- 
istics of the theory describing the system under consideration, while the opposite 
case means a difficulty. 

The questions of passivity is met in the Gaussian white noise problem |lUj and in 
the theory of interacting harmonic oscillators |11[I12) . In this work we demonstrate, 
using quantum transport equations, that the problem of a non-passive reservoir 
occurs in quantum dots. We shall identify a significant source of such non-passivity 
with certain terms in the electron-phonon interaction operator of the system. These 
terms are those, which play a leading role in the creation of the polaron well of the 
charge carrier in a dot. The existence of the overheating in quantum dots shows 
that the formulation of the phonon kinetic equation need to be approached with 
care. We shall show a possibility of a partial elimination of this effect. This will 
be done by changing the basis set of electron and phonon cigcnstates, switching 
from the basis from noninteracting electron and phonons to a set of dressed states 
of electron and phonons, which assumes that the crystal lattice accommodates to 
the presence of electrons in quantum dot. The transition to the new dressed basis 
is achieved by performing the well-known canonical Lang-Firsov transformation 
1181 of the Hamiltonian of the system, which exactly diagonalizes a certain part of 
the full electron-phonon Hamiltonian, including those terms of the electron-phonon 
coupling, which lead to a formation of the electronic polaron well. The result of 
the transformation then is a Hamiltonian for new quasiparticles interacting by a 
new interaction operator. The advantage of this new picture then is that the new 
quasiparticles do not tend to develop quickly a polaron potential well around them. 
The phonon kinetic equation may then become free of corresponding spurious effects 
like an excess heating of phonons. 

Considering electronic relaxation and phonon generation in the system of a 
single quantum dot with a single spin-less electron and with the influence of holes 
in the valence band states of quantum dot completely neglected 0111111, we shall 
demonstrate the existence of the phonon overheating effect in quantum dots and 
show how to eliminate it to a considerable degree. We also show that the electronic 
relaxation characteristics, given by the electronic transport equation, remain nearly 
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unchanged by the treatment of the overheating effect with help of the Lang-Firsov 
transformation. 

2 Hamiltonian 

The electronic eigenstates of this single electron in a quantum dot will be ap- 
proximated by electronic eigenstates in a cubic shaped quantum dot, with infinitely 
deep potential. The present model assumes only two nondegenerate electron states 
il^niy) in this dot, namely the ground state {n—Q) and one of the lowest-energy 
triple-degenerate exited states (ri=l). The bath of LO phonons will be approxi- 
mated by the bulk modes of LO phonons of the full sample. The Hamiltonian of 
the whole system is: 

H = He + Hph + Hi, (1) 

in which ~ Si=o i ^^''i''''*' '^i^h Q being annihilation operator of electron in 
state i. We set the energy of the ground state Eq = 0. Hph = J2q ELob^bq, with 
6q being LO phonon annihilation operator in state with phonon momentum q. The 
interaction of electron with modes q of phonons, with phonon energy E^Oy is 

Hi^ Ag<^{n,m,q){bq~b\)c+Cm. (2) 

q, m,n— 0,1 

The Frohlich's coupling Hi between the electron and LO phonons contains 
the coupling constant Aq, which is Aq — {—ie / q)[ELo{i'^Zc ~ '^vi^)Y^'^{'^^o^)^^^'^ ^ 
where Kqo and kq are, respectively, high-frequency and static dielectric constants, 
Eq is permittivity of free space, — e is the electronic charge, q = | q |, and V is volume 
of the sample. $ is the form-factor, <i>(n, m, q) = J (i'^rV'i^(r)e*'''"V'm(r)- 

This Hamiltonian consists of the electronic subsystem, iJg, the bath Hp^ and 
their coupling Hi. A remarkable feature of this Hamiltonian is the presence of 
two different components of interaction operator. This operator can be written as 
Hi = h'^^ -\- Hf\ The transverse part h'^^ contains that part of iJi, in which 
the electronic orbital indexes n, m are equal to each other, which meaning is that 
electron emits or absorbs phonon without changing orbital state. The longitudinal 
term h'^^ contains terms, in which the electronic orbital indexes are different, 
expressing transition between different orbital states, during emission or absorption 
of a phonon. The former term H^"^ plays a role in processes of development of a 
polaron well, after the electron is placed into an orbital. In the formulation of 
phonon quantum transport equation for the distribution function < 6q6q' > wc 
do not take into account the process of the development of such a polaron well in 
the form of a deformation of the potential in which the harmonic oscillators of the 
phonons move in the presence of the charge carrier. The deformation of the phonon 
field could be described by a quantity like the average < 6q > , for which we are not 
currently considering an appropriate transport equation. Therefore, it would not 
be unexpected if the term Hf' led to a conflicting effect in the standard phonon 
transport equation for the density matrix < 6qfoq' >. The transverse term of the 
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electron-phonon coupling may contribute to this effect already in the first order of 
perturbation, while the longitudinal term i/^''' can contribute to such an effect only 
at higher orders. 

We need to describe the electron-phonon system with having the phonon bath 
passive. This means that energy does not permanently come into or leave the bath 
in the long-time limit. The bath energy is related to the total phonon number 
given by the statistical distribution < b^bq' >. Therefore, it may be meaningful to 
change the basis of the electron and phonon states to such a basis, in which the 
lattice deformation is already included and the lattice does not display a tendency 
to build it up any more, being already accommodated to the presence of electron. 
This representation may be achieved upon performing the well-known Lang-Firsov 
canonical transformation, which is known to provide exact diagonalization of the 
independent boson model [Hi, which in our case corresponds to the Hamiltonian 
Hib=H-H['\ 

We shall show the result of the Lang-Firsov transformation of the Hamiltonian 
H and we shall also present the transformation of the remaining longitudinal part 
of the coupling h[''\ We shall calculate the magnitude of the permanent phonon 
generation in the phonon bath within the quantum transport equation based on 
the full interaction operator Hi. Then we shall show how much the Lang-Firsov 
transformation helps to suppress the phonon overheating effect. 

With help of canonical Lang-Firsov (LF) transformation S we introduce new 
particle operators of electron. An, and of phonons, Bq, by relations c„ — S'A„S'+ 
and 6q = SBqS+. The operator S = , and = - Y.j,cS-^qHh 3, (\)/ELo){Bq + 
B^q)A^ Aj. This transformation diagonalizes exactly the independent boson part 
of the full Hamiltonian, Hjb- Therefore, the full Hamiltonian has now the form 

H = Y^iEn - ^)AUn + E ELoB+Bq - ^AtAiA+Ao + (3) 

where kqi = ^X^q I -^q P -Re[<i>(0, 0, q)<i>(l, l,q)], while the electronic polaron well 

depth is ann/ Elo, with a„„ = X)q I ?\ q) P- ^i'"* is the longitudi- 

nal part of the electron-phonon interaction transformed by the canonical transfor- 
mation S. This operator comes out from the transformation as i/}'"* = Vi + V2. 
Before specifying these two operators, let us introduce an abbreviation, 7j,q — 
Aq^{j,j, q)/ELo, which is a quantity of the first order in the electron-phonon cou- 
phng Aq. The operator V2 reads: 

= 2 5^A,$(0,l,q)7,,-qA+Aie'3 + /i.c., (4) 

where (3 = X]p(7i,p ~ 7o,p)(-Sp + B^^). Let us note that this part of the electron- 
phonon interaction is at lest of the second order in the coupling constant Aq. In 
the lowest order in the electron-phonon coupling constant it gives a phonon-less 
electronic level mixing. We shall neglect this effect. Because in V2 the electron- 
phonon coupling starts in the third order of the expansion in powers of the coupling. 
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we shall neglect V2 in comparison with Vi, which reads: 

Vi^Yl 1' ~ B+^)^oAie^ + h.c. (5) 

q 

The effect of the Lang-Firsov transformation has thus been concentrated in the 
latter formula into the exponential factor, which now contains a complicated de- 
pendence on the phonon operators. The advantage of making the LF transformation 
is also in that we can further simplify the electron-phonon operator Vi . The expan- 
sion of Vi in powers of the coupling starts with the first power of Aq . The lowest 
order term corresponds to putting (3 equal zero. Confining ourselves to keeping only 
the leading order terms in Vi, we estimate the influence of the exponential factor. 
Namely, we approximate e^^ by an average of this quantity < >, calculated 
with canonical ensemble at a temperature T. Calculating this average fXBI we get 
< el^ >= exp[-^p I 7i,p - 7o,p P (< B+Bp > +1/2)]. Here < B+Bp"> has the 
meaning of an average number of phonons with wave vector p. We shall first of 
all assume that this number of phonons is small, comparable with the factor 1/2. 
For the second, we see that the order of the square of difference of the coupling 
constants 7„,p is given approximately by the polaron constant ^Sli which is about 
0.07 in GaAs. From these reasons we approximate the factors by the value of 1. 
After having made these approximations we can see that the longitudinal part of 
the coupling has not changed upon making the LF transformation. 
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Fig. 1. Long-time limit (steady state) of occupation TVo 
of electronic lower-energy state at given temperature Tlo of LO-phonons. Dashed 
line is calculated without LF transformation. 



We also neglect the polaron well depth component of new electron energy, 
ann/ELOi in ^ and the term proportional to kqi in the same equation. These 
terms are small on the scale of the electron energy level separations, and the optical 
phonon energy, in cases we are interested in. With the above made approximations. 
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the resulting Hamiltonian of the system, after making LF transformation, is the 
operator identical with the original Hamiltonian H of equation with the only 
change that the transverse coupling has been dropped. In other words, when 
estimating numerically the electronic relaxation in the dressed electron and phonon 
basis, the original Hamiltonian |^ can be used, with the transverse terms put equal 
zero. The removal of the transverse terms from the original Hamiltonian is done 
simply by putting the form-factors $(n, n, q) equal zero for n = 0, 1. 

3 Electron relaxation and phonon heat generation 

In order to demonstrate the overheating effect, we calculate first the heat gener- 
ation in the phonon system with using the full original Hamiltonian ^ . The kinetic 
equations for the electron was derived earlier 12 |H] with help of either noncquilib- 
rium statistical operator or nonequilibrium Green's functions 6 , 5 . The derivation 
of the kinetic equation for the phonon system can be obtained in the same way. 
The state of the phonon system should be generally described with help of the 
density matrix < b^bqi >, or equivalently with help of the Wigner function. This 
generality would be good for expressing correctly the generation of phonons in the 
area of the quantum dot only. In order to reduce the computing demands as much 
as possible, while conserving reasonably much the prcdictivity of the model, we 
reduce somewhat the complexity of the system of phonon modes. In that sense, 
in order to get rid of the off-diagonal elements of the phonon density matrix, we 
identify the volume of the basic area of the sample with the volume of the quan- 
tum dot. In addition to this, we formally reduce the number of the modes, with 
which the electron interacts, to just a single one, with phonon annihilation operator 
b, the interaction constant of this mode with the electron being chosen such that 
the coupling of electron to this single mode provides the same relaxation rapidity 
as the original full set of LO phonon modes. We shall publish the argumentation 
leading to this simplification elsewhere. The kinetic equations for the electron dis- 
tribution function in the excited state, < c^Ci >, is derived in the self-consistent 
Born Approximation to electronic self-energy and in the random-phase approxi- 
mation to the phonon self-energy. In the right hand sides of the kinetic equations, 
and in the equation for the self-energy, the phonon propagator is approximated 
by undressed phonon propagator. Instantaneous collisions approximation is used. 
Standard Kadanoff-Baym ansatz is used, which step is regarded as suitable enough 
for expressing the long-time limits of development of the system's statistics. Other 
details of the derivation of the kinetic equations will also be described elsewhere. 

We shall first show the electronic relaxation at a given temperature of the LO- 
phonon reservoir. We use the material constants of GaAs, with energy of LO phonon 
of 36.2 meV. The lateral size of the cubic dot is taken to be 16.4 nm. This size 
corresponds to the electronic energy level separation of 62.7 meV. The electronic 
transport equation calculated with the help of the full Hamiltonian in the original 
form 1^ was solved first to find the electronic distribution among the two levels at 
a given lattice temperature. In Fig.^the dashed line shows the lattice temperature 
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Fig. 2. Plot of generation rate d < b'^b > /dt of LO-phonons as a function of temperature 
Tlo of the lattice. The steady-state occupation A'^o of the lower-energy electronic state is 

shown in Fig.0 



dependence of the occupation A^o of the lower energy electron orbital, at which the 
relaxation rate dNo/dt is zero. In other words, at this Nq the electronic subsystem 
reaches the long-time limiting steady state distribution. Let us remark that we have 
earlier shown that these values of Nq are well below the values given by the 
Fermi-Dirac distribution of electron in the present electronic two-level system. The 
dashed curve at Fig. |21 shows the generation rate of LO phonons ^^^^^^ at given 
phonon temperature and at the steady state electronic distribution given by the 
dashed line in Fig.^ One can see that at these conditions the phonon system does 
not have a steady state. It permanently heats up with the rather high speed of about 
1.5 to 1.8 phonons per picosecond (the dashed hne). We ascribe this fast heating 
to the role of the transverse coupling terms in the original Hamiltonian , as dis- 
cussed in an earlier section. From the above given reasons we go over to a different 
representation, using dressed electron and phonon basis, writing quantum trans- 
port equations for electron and the phonon transformed with help of Lang-Firsov 
transformation. In this new representation we use an approximate Hamiltonian, 
which equals the original operator Q with the transverse terms dropped. When 
the same characteristics of electron and the phonon are calculated, we get the full 
lines in both Figs. Q and |21 We see that in the new basis the phonon heating is 
considerably reduced, namely by the factor of about 7, with the present simple tool 
based on LF transformation. 

Summing up, let us note the remarkable numerical result for the long-time limit 
of electronic distribution at given lattice temperature: this long-time limit prac- 
tically does not depend on whether the calculation is done with or without the 
Lang-Firsov transformation procedure applied. This helps us to make a conclusion: 
the use of the original representation in undressed states, using the full Hamiltonian 
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l[T]l having the property of the phonon overheating effect, appears to give us plausi- 
ble data about the electronic relaxation in the quantum dot. From this reason, the 
earlier made conclusions about the electron relaxation effects |Hlin]) may remain 
valid. As for the phonon kinetic equation, the transformation to the dressed basis 
with LF transformation helps to reduce the overheating effect considerably. 
Acknowledgement. This work was supported by institutional project AVOZ10100520. 

References 

[1] A. D. Yoffe, Adv. Phys. 50 (2001) 1-208. 

[2] K. Krai and Z. Khas, phys. stat. sol. (b) 208 (1998) R5-R6. 

[3] H. Tsuchiya and T. Miyoshi, J. Appl. Phys. 83 (1998) 2574-2585. 

[4] O. Verzelen, G. Bastard and R. Ferreira, Phys. rev. B 66 (2002) 081308. 

[5] E. M. Lifshitz and L. P. Pitaevskii, Physical kinetics, Butterworth-Heinemann, 1981. 

[6] D. N. Zubarev, Neravnovesnaya statisticheskaya termodinamika, Nauka, Moskva 1971 
(in Russian). 

[7] K. Krai and Z. Khas, Phys. Rev. B 57 (1998) R2061-R2064. 

[8] K. Krai, P. Zdenek, Z. Khas, Nanotechnology, IEEE Transactions on, 3 (2004) 17-25. 

[9] K. Krai, P. Zdenek, Z. Khas, Surf. Sci. 566-568 (2004) 321-326. 

[10] H. Haug and A.-P. Jauho, Quantum Kinetics in Transport and Optics of Semicon- 
ductors, Springer, Berlin 1998. 

[11] Th. M. Nieuwenhuizen and A. E. Allahverdyan, Phys. Rev. E 66 (2002) 036102. 

[12] A. O. Caldeira and A. J. Leggett, Ann. Phys. 149 (1983) 374-456. 

[13] G. D. Mahan, Many-Parttcle Physics, 2nd ed.. Plenum Press, New York 1993. 

[14] J. Callaway, Quantum theory of the solid state, Part B, Academic Press, New York 
1974. 

[15] K. Krai and P. Zdenek, Physica E, to be publ. 



8 



A 



Czech. J. Phys. 56 (2006) 



